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Abstract 

Recent discovery of mesoscopic electronic structures, in particular the carbon nanotubes, made 
necessary an investigation of what effect may helical symmetry of the conductor (metal or semi- 
conductor) have on the persistent current oscillations. We investigate persistent currents in helical 
structures which are non-decaying in time, not requiring a voltage bias, dissipationless stationary 
flow of electrons in a normal-metallic or semiconducting cylinder or circular wire of mesoscopic 
dimension. In the presence of magnetic flux along the toroidal structure, helical symmetry couples 
circular and longitudinal currents to each other. Our calculations suggest that circular persistent 
currents in these structures have two components with periods 3>o an d &o/s (s is an integer spe- 
cific to any geometry). However, resultant circular persistent current oscillations have 3>o period. 
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Aharonov and Bohm showed that, contrary to the conclusion of classical electrodynam- 
ics, there exists effects of the potentials on the charged particles even in the region where 
all fields vanish. This effect has quantum mechanical origin because it comes from the in- 
terference phenomenon. The well-known manifestation of the Ahoronov-Bohm (AB) effect 
is the oscillation of electrical resistance and the periodic persistent currents in the normal 
metal rings and mesoscopic rings threaded by a magnetic flux. This current arises due to 
the boundary conditions imposed by the doubly connected nature of the loop. Therefore, 
electronic wave function and then any physical property of the ring is a periodic function of 
the magnetic flux with a fundamental period $o- in particular, flux dependence of the free 
energy implies the existence of a thermodynamics (persistent) current. 

Persistent currents in mesoscopic rings: Persistent currents in mesoscopic systems was 
first predicted by one of the authors^ and later discovered by Buttiker 2 et al. A number 
of key experiments also confirmed the existence of persistent currents in isolated rings^. 
In the presence of magnetic flux ($) applied at the center of the ring, we consider a one 
dimensional ring of circumference L r = 1-nr = NA. Here N is the number of lattice points 
and A is the lattice spacing. Tight-binding Hamiltonian reads 



where to is the hopping amplitude between the nearest-neighbour sites for an undistorted lat- 
tice and operators a n (a+) annihilates (creates) an electron at site n. a is the corresponding 
phase change which can be expressed in terms of AB flux ($) 



where $o = hc/e ~ 4.1 x 10 -7 G cm 2 is the flux quantum. 

Corresponding eigenvalue spectrum and the persistent currents (variation of free energy 
with the magnetic flux) are both periodic in $ with a period $o- If we ignore spin of the 
electron, ground state energy and the total current flowing along the ring can be written as 



JV 




(1) 



n=l 




(2) 





n 



n 



de n (<&) 




(4) 



2 



where Iq = AcKt^/N^Q is the current amplitude, and the summation is over number of 
electrons, N e , for each value of flux. 

Persistent Currents in Helical Structures: Persistent currents was believed to be a spe- 
cific property of isolated systems for a long time^. However, theoretical studies suggest that 
persistent currents should also exist in connected rings- In the presence of both longitudinal 
and transverse flux, existence of transverse persistent currents in doubly connected meso- 
scopic rings is known^. Moreover, transverse currents may contribute to an experimental 
observation of longitudinal persistent current and it can substantially increase the amplitude 
of the AB oscillations^. 

AB effect is also shown to be present in toroidal systems^. In this paper, we consider a 
set of identical and connected mesoscopic rings with a circumference of L r = 2irr and each 
having N lattice sites with a lattice spacing of Ai = L r /N. We also assume our helical 
structure has L periods with a periodicity of iV rings as shown in Fig. So, we have a 
toroid of circumference L t = 2ttR and it contains LN rings which are uniformly separated 
by A2 = Lt/LN along the circumference of the toroid. In the presence of magnetic fluxes 
$(a) and $(/?), which are applied through the center of the connected rings and at the 
center of the toroid respectively, we propose two models. 

First Helical Model: In order to have helical symmetry, we allow only nearest-neighbour 
circular hopping between the sites in each ring and vertical hoppings (no cross-hoppings) 
between the nearest-neighbour rings along the toroid respectively. Assuming the tight- 




FIG. 1: Left: Toroid of LN = 3*3 = 9 rings are connected by vertical hoppings and is 
applied at the center. Right: Each ring has N = 3 sites and <5(a) is applied at the center. 

binding model for electron transport, system Hamiltonian can be written as 

L N N 
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where t% is the hopping amplitude between the nearest-neighbour sites in the ring, t 2 is the 
vertical hopping amplitude between the nearest-neighbour rings, t± is the special hopping 
amplitude in vertical direction and £3 is the special circular hopping amplitude which con- 
nects two special vertical hoppings as shown in Fig. ^ Note that in order to study helical 
symmetry with this Hamiltonian, it is necessary to have t 3 ^> t\ and £4 ^> t 2 . Operator 
a ikn ( a tkn) annihilates (creates) an electron at period 1, ring k and site n. n N represents 
n(modN), 8^ is the Kronecker delta, h.c. is hermitian conjugate and a and 3 are the corre- 
sponding phase changes between nearest sites in a particular ring and between nearest rings 
along the toroid respectively. They can be expressed in terms of AB flux (<&) as 

$(a) „ MB) 
iV$ LiV$ 

Eigenf unctions of the hamiltonian, H = J2i kn^i'k'n' hi kn ,i'k'n'( l t kn a i'k'n' can be written 
as an expansion of ip = J2i,k,n^ikn a ik n \^ > where |0 > denotes the vacuum state. Expansion 
coefficicents C ikn satisfies Ewv hikn,Vk'n'C Vk 'n' = E lkn Ci kn . According to Fig. □ system 
exactly repeats itself after translation along the toroid by N rings. Therefore Bloch theorem 
applies and it gives Ci + i tk ^ n = e^Ci^n where 7 = with s — 0, 1, L — 1. 

Bloch theorem partly digonalizes matrix Hi kn ^i k i n i by quantized values of 7 with a reduced 
Hamiltonian matrix elements given by 
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Matrix H kn k > n > should be diagonalized numerically which suggests that instead of k, n, 
we introduce r = N(k — 1) +n with r changing from 1 to N 2 . For a given value of a and 
3, Hamiltonian has total of LN 2 eigenvalues since Hk nj k>n'['l) should be diagonalized for 
each value of 7 = ^-s, s = 0, 1, ...,L — 1. Finding these eigenvalues accomplishes unitary 
transformation of creation operators from af kn to a+, which ensures that states ip rr i = 



ai./ 1 > are orthogonal to each other and (i^ r i are canonical Fermi operators, i.e. [a+ , a r v]+ = 
o SS fS rr ' . 

For a given number of electrons in toroid, N e , we calculate minimal energy which is sum 
of lowest N e out of LN 2 eigenvalues and also calculate total persistent currents both along 
the toroid and along the rings by variation of free energy of the system with the magnetic 
flux. 




FIG. 2: Contour plot of I c (a) vs. Flux in the first helical model. Parameters are t± = = 10, 
t 2 = t 4 = 1, JV = L = 10 and N e = 500 (half filling). 
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FIG. 3: Contour plot of Ii{(3) vs. Flux in the first helical model. Parameters are t\ = t% = 1, 
t 2 = t 4 = 100, N = L = 10 and N e = 500 (half filling). 



Total persistent currents along the rings (circular) and the toroid (longitudinal) are per- 
pendicular to each other. In Fig. |21 we choose hopping parameters such that probability of 
finding electrons inside the rings is much more than along the toroid. In this limit, I c (a) 
dominates Ii(j3) — > and the coupling between these currents is very small (negligible) as 
shown in the figure. In the opposite limit, where electrons have more probability of being 
along the toroid than inside the rings, Ii{f3) dominates Ii(a) — > and the coupling is also 
very small as shown in Fig. El 

In order to understand the mixing of both symmetries, we choose another set of param- 



5 



eters such that probability of finding electrons along the helical path (Fig. ^) is much more 
than finding it elsewhere in the toroid. This case is shown in Fig. 0] and |S| As opposed to 
the previous limits, we find that mixing symmetries has cross-effect on the system and both 
currents coupled to each other. However, both currents in different directions are periodic 
in $ with a period of $o as expected in all cases. 




FIG. 4: Contour plot of I c (a) vs. Flux in the first helical model. Parameters are t\ = t<i = 1, 
t 3 = t 4 = 10, N = L = 10 and N e = 500 (half filling). 
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FIG. 5: Contour plot of Ii((3) vs. Flux in the first helical model. Parameters are t± = t2 = 1, 
t 3 = t A = 10, N = L = 10 and iV e = 500 (half filling). 



Second Helical Model: In order to satisfy necessary boundary conditions for the geometry 
of the toroid ring m = LN + 1 should coincide with ring m — 1. As shown in Fig. 
we fix the position of the first ring and rotate rest (LN — 1) of them by 2irs/LN where 
s = 0,1, ...,LN — 1. Note that each value of parameter s specifies different geometry. In 
this model, we considered all possible circular and cross hoppings both inside and between 
the rings. 

In the tight-binding approximation, system Hamiltonian becomes 

LN N 



m=l n=l 
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FIG. 6: Left: Connected rings in the first model with <£(a) applied at the center. Right: Rotated 
rings by angle 2tts/LN for s = 1, N = 3 and L = 3. 
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where *i is the hopping amplitude between the nearest-neighbour sites in a ring and *2 is 
the cross-hopping amplitude between the rings. Operators a mn (a+ n ) annihilates (creates) 
an electron at ring m and site n. a and (3 are the corresponding phase changes between sites 
in a particular ring and different rings in the toroid respectively. They are given by Eq. |3 
Hamiltonian can be diagonalized by discrete Fourier transformation from a mn in site to 
b q k in momentum representation as 



1 U J(knA 1 +qmA 2 ) 



(9) 



where k = and n = 0, 1, 2, N — 1 and 

diagonal form, Hamiltonian becomes 
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q = L ^' A m and m = 0, 1, 2, LN — 1. In 
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Eigenvalues of this Hamiltonian are periodic in <£> with a period of $0 an d they are given by 
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Corresponding total persistent currents along the rings and the toroid (circular a and 
longitudinal f3 currents) which are periodic in $ with a period of $0 are perpendicular to 
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each other and they are given by J c ($(a)) = J2m,n ^frffi where summation is over number 
of electrons, N e . Ignoring the spin of electrons, for each value of flux we have 
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where I\ = — 4c7rt 1 /A^$ an d I2 — — 4c7it2s/ LNQ are the current amplitudes. 

Discussion and Conclusion: The geometric structure determines the electronic structure 
and thus the characteristics of the persistent current oscillations. In this paper, we study 
symmetry mixing and cross-effects in a toroidal system which is threaded by magnetic flux 
both along (a) and inside (J3) the structure. We consider a set of connected-mesoscopic 
rings and model helical symmetry by restricting hopping directions in our models. The 
electronic structure calculated from the tight-binding model is given in Eq. ^2 for second 
helical model, however, we can not solve it for the first helical model and instead evaluate 
it numerically. Since magnetic flux $(a) and $(/3) are in perpendicular directions, circular 
and longitudinal currents also flows in perpendicular directions. 

In both models we consider electron transport inside the rings, however, we propose 
two models with hoppings in different directions between the rings. In the first model, we 
allow only vertical hoppings between the neighbouring rings. Since the system is periodic 
along the toroid, we include the Bloch condition and solve the final Hamiltonian matrix 
for energy eigenvalues and persistent currents numerically. Our results show that mixing 
perpendicular magnetix fluxes couples perpendicular currents with each other and both 
circular and longitudinal currents are periodic in $ with a period $o- 

In the second model, we consider electron transport with cross hoppings between the rings. 
This coupling between perpendicular a and (3 magnetic fluxes yields an extra component 
to the total circular current (fT2"j) with a period $o/s. Since s is a positive integer, total 
circular persistent currents have period $0 as expected. In the special case, for s = 0, all 
cross-hoppings are indeed now vertical hoppings and we recover the result of first model 
together with circular currents which have period $o- We also note that Eq. ^Jand Eq. EI] 
are in agreement with Eq. 0]in the limits when N — > l,s — and t 2 — > 0. 
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Extra component (coupling) of circular persistent currents appears in both models. These 
currents are vanishingly small in the limit of large number of rings, L — > oo, as expected. 
Note that, extra circular current component is due only to $(a) and the presence of $(/3) 
results only in longitudinal persistent currents along the toroid. Our both model results 
are also in agreement with Lin- et al.. They showed that perpendicular $(/3) through the 
carbon nanotube toroidal structures results in persistent current oscillations with a period 
$o- To conclude, our calculations suggest that circular persistent currents in structures with 
helical symmetry have two components with periods $o an d $o/ s - However, total circular 
persistent current oscillations have <3>o period. 
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